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Big Bang nucleosynthesis in the presence of decaying relic neutral particles is examined in detail. 
All non-thermal processes important for the determination of light-element abundance yields of ^H, 
■^H, "^He, ^He, ®Li, and ^Li are coupled to the thermonuclear fusion reactions to obtain comparatively 
accurate results. Predicted light-element yields are compared to observationally inferred limits on 
primordial light-element abundances to infer constraints on the abundances and properties of relic 
decaying particles with decay times in the interval 0.01 sec < tx ^ lO^^sec. Decaying particles are 
typically constrained at early times by *He or ^H, at intermediate times by ®Li, and at large times 
j-S^ , by the '^He/^H ratio. Constraints are shown for a large number of hadronic branching ratios and 

^— ^ ■ decaying particle masses and may be applied to constrain the evolution of the early Universe. 
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I. INTRODUCTION 



The epoch of Big Bang nucleosynthesis (BBN) is one of the furthest reaching back probe of early cosmological 

conditions. It has thus been invaluable, for example, for the realization that baryonic matter contributes only a small 

fraction to the present critical density, thus providing independent evidence for the existence of dark matter, or that 

extensions of the standard model of particle physics which include further light degrees of freedom are observationally 

l/~) ■ disfavored. Due to its usefulness BBN has thus been analysed in many variants (inhomogeneous BBN, BBN with 

O^ ' leptonic chemical potentials, BBN with antimatter domains, BBN with decaying particles or evaporating primordial 

">! i black holes, BBN with Brans-Dicke gravity, etc. - for reviews the reader is referred to Ref. [l|). These studies have 

rj^ [ been important for the realization that the Universe must have been fairly close during the BBN epoch (^^ 1-lO^sec) 

f— ^ . to that described by the standard model of BBN (SBBN). An SBBN scenario assumes a homogeneous baryonic 

"^ I plasma free of defects or decaying relics and with close-to-vanishing leptonic chemical potentials at constant baryon- 

f^ ■ to-entropy. With the recent accurate determination of the fractional contribution of baryons to the critical density 

I ' Vtbh"^ = 0.0223+00009 {h is the Hubble constant in units 100 km s^^Mpc^^) by the WMAP sateUite mission [I 

j~^. predicted primordial light-elements in a SBBN scenario may be compared with those observationally inferred. This 

f— I I comparison is good when ^H/H is considered, fair when ^He, and inconclusive when ^He is considered, but factor 

• • . 2 — 3 discrepant in the case of ^Li. It is not clear if this latter discrepancy is due to stellar ^Li destruction effects, 

.5^ ■ either in Pop II [3|, Ij, la, l6| stars where it is observed, or in a prior Pop III generation 7], or indeed due to a required 

\^ I modification of BBN, such as a population of decaying relics [8|, |9|. This discrepancy, however, is not subject of the 

^ ' present paper. 

The current study presents in detail analysis and results of BBN with decaying neutral particles. Such scenarios 
may easily emerge, for example, in supersymmetric extensions of the standard model, where relic gravitinos produced 
at high temperatures decay during BBN [lOl. Similarly, when the gravitino is the lightest supersymmetric particle, 
next-to-lightest supersymmetric binos or sneutrinos may decay into gravitinos duringor after BBN. Many such BBN 



studies have been presented before, from the pioneering works in the eighties llll . Il2l Il3l . 114 Il5l | to improved and 
updated analysis in the nineties [l^, lul, Il8| up to very recent publications [3, llSi l20!, l2l| . Nevertheless, except for 
Ref. 8, 20, 21] none of these works were able to properly treat the hadronic decay of a relic particle during BBN. 

The purpose of this paper is twofold. It is intended to present a catalog of constraints on hadronically and 
electromagnetically decaying neutral particles over a wide range of decay times (10~^-10^^sec) as well as hadronic 
branching ratios and relic particle masses. Prior studies have presented results only for a very small number of 
hadronic branching ratios. It should be noted that due to nonlinearities, either at smaller decay times (r < lO'^sec) 
when thermal nuclear reactions still proceed or due to interference of electromagnetic destruction and hadronic 
production, for example, a simple interpolation between results at different hadronic branching ratios often leads 
to erroneous or inaccurate results. The second purpose of this study is to present the detailed analysis underlying 
the results presented in Ref. [3, |9|, 122, l23| . In general the philosophy in the analysis is to stay as close as possible 
to experimental data. Furthermore, a large number (93, cf. Table HJ of processes are included, all those which the 
author believes to play a role at the more than a few per cent level. Particular attention is given to an as accurate as 



possible determination of the synthesized ^Li in such scenarios. This is also of importance as ^Li has been recently 
observed in surprisingly large abundance in about 10 Pop II stars [24| and is often a byproduct of decaying [IJ, |l7| 
or annihilating 2M ^^^^c particles. Finally, it is noted that the present limits do not apply to decaying electrically 



charged relic particles with tx ^ 100 sec as in this case, the existence of bound states between nuclei and the relics 



may significantly change nuclear reaction rates 25 1 



The outline of the paper is as follows. In Sec. II a brief overview over the physics of BBN in the presence 
of electromagnetically and hadronically decaying particles is presented. Sec. Ill discusses observationally inferred 
primordial light-element abundances and derives the constraints on those as applied in the analysis. Sec. IV then 
presents the constraints which apply to the abundance of putative relic decaying particles, whereas conclusions are 
drawn in Sec. V. For the in details interested reader and for the benefit of the author, a number of appendices 
summarizes the analysis. This includes general considerations (App. A), the numerical procedure (App. B), employed 
kinematic relations (App. C), non-thermal electromagnetic interactions (App. D), non-thermal hadronic interactions 
(App. E), and the thermal nuclear reactions (App. F). 

II. OVERVIEW OF THE PHYSICS OF BBN WITH DECAYING PARTICLES 

This section gives a brief and schematic account of the modifications to BBN when energetic particles are injected 
due to, for example, the decay of relic particles. The decay channel of the relic determines the energy and nature of 
the injected primaries. A BBN calculation with decaying particles requires the detailed study of the thermalization of 
these primaries in the primordial plasma, as well as that of secondaries produced during the thermalization process. 
Since for one energetic primary several tens/hundreds of secondaries may be produced such calculations are sometimes 
also referred to as cascade nucleosynthesis calculations. In addition, thermal nuclear reactions operative at early times 
(r < lO'^sec) have to be followed. One distinguishes between primaries being either hadronic (nucleons, antinucleons, 
mesons), electromagnetic (photons, e ), or inert (neutrinos or other weakly interacting particles). This distinction 
is useful as the thermalization of hadronically and electromagnetically interacting particles is very different and thus 
impacts the light-element yields in different ways. Both have nevertheless in common that over the whole range of 
decay times considered here (lO^^-lO^^sec) thermalization occurs very rapidly on the Hubble scale. Redshifting of 
energetic particles has therefore not to be considered. In contrast, inert particles, such as neutrinos have typically no 
further interactions with the plasma and thus also do not effect BBN yields. However, even neutrinos, if energetic 
enough {^ 100 GeV) and injected very early on (r ~ 1 sec) may not be regarded as inert (cf. App. E.l), but since 
their effects are typically subdominant compared to hadronically and electromagnetically interacting particles, they 
are not considered in the present study. 

Hadronically interacting primaries may affect light-element yields in the whole above quoted decay time range. 
Charged mesons have an effect only during early times (r ~ 10~^ — lO^sec) due to charge exchange reactions [15| 
(Rec. 34-43 Table U App. E.l) with nucleons, converting mostly protons to neutrons and thereby increasing the 
"'He abundance. At larger decay times their effects are negligible since they mostly decay before reacting with 
nucleons. At shorter decay times thermal weak interactions very quickly re-establish the SBBN ncutron-to-proton 
ratio. Antinucleons may have an effect at all times, though their effect is the most pronounced during the time 
interval r ^ 10~^ — lO^sec (App. E.2). As they are more likely to annihilate on protons they also tend to increase the 
neutron-to-proton ratio. At times r ^ 200 sec a significant fraction Yp w 0.25 of all baryons are bound into helium 
nuclei. Annihilation of antinucleons on helium nuclei may leave ^H,^H, and ■^He as secondaries [11]. However, this 
effect is subdominant as compared to the effects by injected nucleons. 

Neutrons at r ^ 200 sec and protons at r > lO'^sec thermalize to a substantial degree via nucleon-nucleon collisions 
(Rec. 48-73 Table HI App.E.3,E.4) and nuclear spallation reactions (Rec. 74-83 Table HI App.E.5). This is because 
the competing processes for thermalization, magnetic moment scattering of neutrons off plasma e (Rec. 3, App. 
D.l.c), Coulomb stopping of protons by plasma e (Reel, App. D.l.a), and Thomson scattering of protons on 
CMBR photons (Rec. 2, App. D.l.b) loose their efficiency. This holds true except for very energetic nucleons as 
well as protons in the tens of MeV range and is due to ever decreasing numbers of e^ and energy of CMBR photons. 
These results may be observed in Fig. [TJ Each nucleon-nucleon scattering produces another energetic nuclcon. The 
injection of an energetic 100 GeV nucleon may therefore lead to the production of several tens of 1 GcV nucleons. 
Their collective impact on the light-element yields via spallation of ''He to produce ^H,'^H, and ■^He, is far greater 
than that of a 100 GeV injected antinucleon '14]. The latter undergoes maximal 2-3 scattering before annihilating, 
with the respective probability of spallation of '^He being compartively small. The injection of nucleons is therefore 
always important. At early times this is due to an increase of the neutron density (r ~ 1 — lO'^sec) leading to either an 
increased ''He (r < 200 sec), increased ^H (r ~ 200 — lO'^sec), or decreased ^Li (r ~ lO'^sec) abundance. Somewhat 
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FIG. 1: The quantity _E(Zjvid_B/d2;jc)~^ as a function of nucleon energy E, where /jv is the nucleon mean free path and Id^E/dxIc 
is the nucleon energy loss per unit path length due to 'continuous' energy losses such as multiple Coulomb scattering. When 
-B(Zjv|d_E/da;|c)~^ ^ 1 nucleons predominantly loose energy due to nucleon-nucleon scattering and spallation processes, whereas 
in the opposite case nucleons loose their energy predominantly continuously via multiple electromagnetic scatterings on e* and 
photons (cf. Eq. B.l). Only in the former limit nuclear cascades occur. The quantity is shown for neutrons (green - dashed) at 
temperatures T = 90 and 30keV and protons (red - solid) at temperatures T = 30, 10, 1 and 0.01 keV as labeled in the figure. 



0.01 



0.001 



„^ 0.0001 



le-05 



le-06 



lOkeV 
5keV 

IkeV 
O.lkeV 




10 



100 



E (MeV) 



FIG. 2: Probability PeL; of energetic "^H nuclei (red - solid) and energetic ''He (blue - dotted) of initial energy E to fuse on 
ambient '*He nuclei to form ^Li via the reactions ■'H(a,n)^Li and ^He(a,p)®Li, respectively. The figure shows this probability 
for ''H-nuclei at temperatures T — 5, 10, 5, 1 and 0.1 keV and for ^He-nuclei at temperature T — 0.1 keV, respectively. Survival 
of the freshly formed Li nuclei against thermal nuclear reactions is also taken into account as evident by the comparatively low 
PeLi at T = 15keV. The figure illustrates that Psl; dramatically increases at T ~ 5 — lOkeV due to a decrease of the efficiency 
of Coulomb stopping in this narrow temperature interval (see text). 
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FIG. 3; Yields of (from top to bottom) neutrons (red), deuterium nuclei (green), tritium nuclei (blue), and "^He-nuclei (magenta) 
as a function of cosmic temperature T per hadronically decaying particle X ~> qq of mass Mx ~ 1 TeV, where q denotes a 
quark. Note that initially after hadronization of the g-g state on average only 1.56 neutrons result. The remainder of the created 
neutrons at lower temperatures T < 90 keV are resulting from the thermalization of the injected neutrons (and protons) due to 
inelastic nucleon-nucleon scattering processes and "^He spallation processes. Similarly, all the ^H,^H, and ^He nuclei are due to 
"^He spallation processes and np nonthermal fusion reactions (for ^H) induced by the thermalization of the injected energetic 
nucleons. The figure does note include the electromagnetic yields (cf. Fig. Q due to photodisintegration which is inevitable 
even for a hadronic decay since approximately 45% of the rest mass energy of the decaying particle is converted into energetic 
7-rays and energetic e after hadronization and pion decays. 



later, and for the whole range of decay times considered here (r > 200 sec) spallation of ''He by energetic nucleons 
and their secondaries yields large amounts of ^H,^H, and '^He. These trends may be seen in Fig. [H As '^He spallation 
processes by neutrons in the time range r ~ 200 — IC^sec are by far more efficient that those by protons, inelastic 
nucleon-nucleon reactions (Rec. 52-73, App.E.4) which preferentially convert protons to neutrons are also important 
for accurate BBN yield predictions. 

The spallation of ^He to form '^H and '^He is also important as it may lead to the synthesis of ^Li [14] (cf. App. E.6). 
The mass three nuclei '^H and '^He emerge energetic from these reactions with typical energies around 10 MeV. They 
may thus participate in the non-thermal fusions reactions 85-86 to form ^Li. The formation of ^Li by those energetic 
3-nuclei, as well as by energetic ^He resulting from elastic nucleon-^'He collisions (Rec. 87-88) is by far less important. 
Energetic mass three nuclei to the largest part loose their energy efficiently via Coulomb stopping. However, even the 
comparatively small fraction ~ 10^^ (cf. Fig. [2]) which fuses may yield an observationally important ^Li abundances. 
The freshly fused ^Li may survive thermal nuclear reactions only above r > lO'^sec. A particularity in the Coulomb 
stopping (cf. App.D.l.a) which renders it very much less efficient for nuclei with velocities below the thermal electron 
velocity, which is the case for ~ 10 MeV 3-nuclei at r ~ lO'^sec, increases the ^Li production efficiency around 
this time by a factor of ten. As explained in App. E.6 ^Li yields due to energetic mass three nuclei produced by 
"^He spallation are somewhat uncertain and may well be some tens of per cents higher. The uncertainty is due to 
an incomplete knowledge of the high-energy tail of produced mass three nuclei which may well include a backward 
scattering peak. The present study is therefore conservative in estimating ^Li yields. 

Electromagnetically interacting particles may impact BBN yields only at comparatively large decay times (r ^ 
lO^sec). This is due to energetic photons and e^ interacting efficiently with CMBR photons (cf. App. D.2 ). 7-rays 
rapidly pair-produce on CMBR photons (Rec. 4) and electrons inverse Compton scatter (Rec. 5). Comparatively long- 
lived against further interactions are only 7-rays of energies below the e pair production threshold. This threshold is 
a.t E Ri 2.2 MeV at r '--^ lO^sec and increases continuously as the CMBR temperature drops reaching E « 19.8 MeV at 
T ~ lO^sec. These latter energies are the thresholds for photodisintegration of ^H and ''He, respectively. Once 7-rays 
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FIG. 4: Yields of (from top to bottom) '^He (blue), "^H (green), ^H (red), and "Li (magenta) nuclei per TeV of electromagnetically 
interacting energy injected (in form of energetic 7-rays and energetic e ) due to photodisintegration reactions (^H, ^H, '^He, 
and ''Li) and fusion reactions (''Li) as a function of cosmic temperature. 

have dropped below the pair production threshold their interactions close to the threshold are dominated by scattering 
of CMBR photons (App. D.2.d), and further below the threshold by Bether-Heitler pair production (App. D.2.a) and 
Compton scattering (App. D.2.c). A substantial though small fraction / ~ 0.01 may photodisintegrate either ^He 
or ^H [I4I. The reader is referred to Fig. [4] and Fig. [5] for typical light-element production- and destruction- factors 
during electromagnetic decays. When '^He photodisinthegration becomes important (r > lO^sec), photodisintegration 
of ^H is a subdominant effect as more ^H is produced in the relatively larger numbers of ^He-nuclei which are 
photodisintegrated. Energetic '^H and '^He produced in large numbers during the "^He photodisintegration may lead 
via non-thermal fusion reactions to ^Li [l7|, similar to the energetic mass three nuclei produced in ''He spallation 
during hadronic decays. To a less important degree ^Li may also be created by direct photodisintegration of ^Li and 
^Be. 

Thus, the effects of injection of electromagnetically interacting particles on BBN yields are due to photodisinte- 
gration and are only operative at comparatively late times (r > lO^sec). A much simplifying aspect of BBN with 
electromagnetically decaying particles is the fact that results almost always only depend on the total amount of elec- 
tromagnetically interacting energy, rather than the energy and nature of the injected primary. This is due to the 
number of cascade photons and electrons being so large that asymptotically the same state is reached, independent 
of the initial state. Finally it is noted that hadronic decays (such as X ^ qq, with q a quark) generically also lead to 
the injection of electromagnetically interacting primaries, as for example due to tt^'s which decay into two photons. 

III. OBSERVATIONAL LIMITS ON THE PRIMORDIAL ABUNDANCES 

The following provides a compilation of the observational limits on primordial abundances adopted in the present 
analysis. For more detail on observations and interpretations, the reader is either referred to the original literature 
or reviews on BBN nucleosynthesis (cf. |26l. l27|). 



^He 



The primordial ^He abundance is inferred from observations of hydrogen- and helium- emission lines in extragalactic 
low-metallicity HH-regions. The two most recent observers determinations yield Yp « 0.2421 ± 0.0021 28[ and 
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FIG. 5: Number of destroyed (from top to bottom) ''He (blue), ^H (red), ^He (green), and ^Li (magenta) (^Be - light blue) 
nuclei per TeV of electromagnetically interacting energy injected into the primordial plasma at temperature T. For the purpose 
of illustration we have taken the '^Li and ^Be abundances equal at ^Li/¥L — ^Be/H ~ 4.34 x 10"^''. In reality it is the sum of 
both isotopes which is synthesized at (^Li + ''Be)/II ^ 4.34 x 10^'° in a SBBN scenario at n^/i^ ^ 0.02233 with ''Be being 
converted to Li by electron capture at cosmic temperatures T ~ 0.1 — 1 keV. 



Yp w 0.239± 0.002 [29|]. Nevertheless, somewhat smaller Yp « 0.2345 or larger Yp w 0.2443 determinations by the same 
authors, but determined from (partially) different data sets of, have also been cited. The estimates are significantly 
below the SBBN Yp prediction at the WMAPIII baryon density. However, the quoted error bars are purely statistical in 
nature and only some of the many possible systematic errors could be estimated. Systematic errors include underlying 
stellar absorption, ionization corrections, temperature fluctuations, and coUisional excitation. Recent re-analysis of 
the original data by an independent group '30l], has led to a significantly revised Yp « 0.249±0.009 estimate, including 
a substantial increase of the error bars. Nevertheless, even this latter study does not include a study of the majority 
of systematic uncertainties, such that its utility is somewhat questionable. Many of these systematic uncertainties 
would lead again to a decrease of the inferred Yp. In a conservative spirit, however, the present analysis utilizes the 
upper limit as quoted above 



Yp < 0.258. 



(1) 



A lower limit on Yp is not considered as usually of no use for deriving limits on relic decaying particles. 



B. ^H 

The most accurate method for an observational determination of the primordial ^H/H-ratio is believed to be 
due to observations of low-metallicity quasar absorption line systems (QALS). Such gaseous systems, on the line 
of sight between a high redshift quasar and the observer, produce absorption lines in the continuous spectrum of 
the quasar at the redshifted wavelength of the Lyman-a ^H and ^H positions, respectively. If the QALS has a 
simple velocity structure, the depths of the absorption troughs may be used to infer the ^H/H-ratio. As observations 
focus on systems at very low metallicity, the assumption is typically that only very little stellar ^H destruction has 
occurred, and that the determined abundance may be directly compared to the predicted BBN ^H/H-ratio. When 
the currently six best determinations [3]| are averaged one obtains a QALS-infcrrcd primordial ^H/H abundance of 
^H/H « 2.4 ± 0.4 X 10"^ ^. This value compares favorably to that predicted by a SBBN scenario at the WMAPIII 
baryonic density. Nevertheless, caution has to be applied as the intrinsic dispersion between individual central values 



of ^H/H determinations is much larger than most of the quoted error bars for individual QALS. Thus, the highest 
individual 2-cr upper limit and the lowest individual 2-cr lower limit of all six systems span a range 

1.2 X 10"^ < 2H/H< 5.3 X 10-^ (2) 

much larger than the above quoted error bar. Furthermore, a systematic trend of decreasing ^H/H with increasing 
column density Nfj of the QALS may be found indicating the existence of possible systematic errors. Attempting 
to remain on the conservative side the above range is adopted in the analysis for an acceptable ^H/H abundance. 
Conservatism may be in order not only because a recent proposal of a very active stellar Pop III population, which 
could potentially reduce the primordial ^H/H-ratio in certain (high Nh) systems without much trace of associated 
iron or silicon production [7|. 

C. ^He/^H 

Precise observational determinations of '^He/H-ratios are only possible within our galaxy. Since our galaxy is 
chemically evolved, it is difficult to make a straightforward connection to the primordial "^He/H abundance. Since 
■^He, in fact, is known to be destroyed in some stars, but produced in others, observed ■^He/H-ratios may also only with 
difficulty be utilized to obtain an upper or lower limit on the primordial ■^He/H. This is in contrast to the ^He/^H- 
ratio. ^H is known to be always destroyed in stars (in fact is burned into "^Hc), whereas ^He is either destroyed or 
produced. The cosmic '^He/^H-ratio may therefore only increase with time. Thus a determination of this ratio in any 
environment, even chemically evolved, may be taken as an upper limit on the primordial "^He/^H-ratio 13^. From 
a combination of solar wind observations and the planetary gas component of meteorites and Jupiter, it is possible 
to infer the ^He/Hw 1.66 ± 0.05 x 10"^ and ^H/Hw 1.94 ± 0.39 x lO'^-ratios at the time of solar system formation 
independently [33|. Employing the 2-cr upper end of the '^He abundance and the 2-a lower end of the ^H abundance, 
when may derive a conservative constraint of 

3He/2H<1.52, (3) 

for the primordial '^Hc/^H ratio. 

D. ^Li 

For a long time it has been known that the ^Li/H-ratio in the atmospheres of low-metallicity PopII stars is essentially 
constant with metallicity (the Spite plateau) with deduced central values on the plateau falling in the range 1.23 x 
lO""'^" < ^Li/H ^ 1.73 X 10"^" [3J|. As this value is traditionally interpreted to be close to the primordial ^Li/H-ratio, 
a factor 2 — 3 discrepancy between the SBBN predicted primordial ''Li/H at the WMAPIII baryon density (usually in 
excess of 4 x 10^^*^ and that inferred from the stars on the Spite plateau is apparent. This situation has neither been 
changed by a large number of re-evaluations of the ''Li/H-abundance of the Spite plateau, nor convincingly shown to 
be solved by stellar depletion of ^Li. As the ^Li isotope is not very important in constraining relic decaying particles 
(rather, it may resolve the problem), the reader is referred elsewhere for discussion of this problem (e.g. [Tl. M. |26J). 
The present analysis adopts only a lower limit on the primordial ^Li/H 

^Li/H > 0.85 X 10"^°, (4) 

chosen by the 95% confidence level lower limit as derived by Ref. [35|. 

E. ''Li 

The isotope of ^Li can be particularly useful in deriving constraints on relic decaying particles [i7|. Recently, the 
number of claimed preliminary detections of ^Li/^'Li-ratios in low-metallicity Pop II stars has multiplied by a large 
factor [24| (cf. to Ref. [36| for the few former detections) falling in the range 0.03 < ^Li/^Li < 0.07 with average 
^Li/^Liw 0.042 and with no star having ^Li/^Li in excess of 0.1. The origin of this ^Li is unknown. Concerning a limit 
on the primordial ^Li-abundance the situation is somewhat problematic. This is particularly so, since the predicted 
SBBN ^Li is in excess of that observed. If this latter fact is after all explained by stellar ^Li destruction, ^Li would 
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FIG. 6: Conservative BBN constraints on the abundance of relic decaying neutral particles as a function of their life time for 
a Mx ~ ITeV particle with hadronic branching ratio Bh = 1. Limits are given on the contribution the decaying particles 
would have made to the present critical density, Q,xh^ (with h the Hubble parameter), if they would have not decayed. For a 
conversion to constraints on, for example Mxnx/n^ the reader is referred to App. A. The colored regions are excluded and 
correspond to the constraints imposed by the observationally inferred upper limit on *He - orange - (Eq. [T|), upper limit on ^H 
- blue - (Eq. [2]), upper limit on '^He/'^H - red - (Eq. [3|, and lower limit on ^Li - light blue - (Eq. |4ll. Conservative constraints 
derived from ®Li/^Li (Eq. |6} are shown by the green region. The region indicated by yellow violates the less conservative 
^Li/^Li (Eq. [5]) constraint but should not be considered ruled out. Rather, this region may be cosmologically interesting as a 
putative source of ®Li in low-metallicity stars by relic decaying particles. 
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FIG. 7: As Fig. E] but for hadronic branching ratio Bh = 3.333 x 10"^ 
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FIG. 8: As Fig. [6] but for hadronic branching ratio Bh = 10 "*. The region excluded by the lower limit on ^H/H (Eq. [2} is 
indicated by the color magenta. 

also be destroyed. In the best case (destroying as little as possible ^Li), both the ^Li and ^Li abundances would be 
destroyed by the same factor, with their ratio, nevertheless, staying constant. In this case one may adopt 



^Li/^Li < 0.1 



(5) 



as upper limit. This would correspond to a primordial ®Li abundance ^Li/H < 4.34 x 10^^^ (given that the SBBN 
prediction is 4.34 x 10^^" in the present analysis, cf. App.F) already much in excess of the 2-cr upper limit of 
^Li/H < 1.47 X 10~^^ in the well studied star HD84937. However, calculations of ^Li and ^Li destruction in stars 
generally predict more ^Li destruction than ^Li destruction, since ^Li is more fragile than ^Li. In Ref. [31] which studied 
Li depletion in rotating stars a correlation between the amount of ^Li and ^Li destruction was obtained, with, for 
example, factor 2 (4) reduction of the ^Li/^Li ratio occurred for factor 1.6 (2.5) reduction of the ^Li abundance. When 
the (fine tuned) models by Ref. [J] are taken, designed to explain a larger amount of ^Li destruction without spoiling 
the flatness of the Spite plateau (when plotted against stellar temperature), a factor 1.6 (2.5) of ^Li destruction implies 
a factor 1.25 (15.8) reduction in the ^Li/^'Li-ratio. Particularly the last destruction factor would render the ^Li isotope 
not anymore as useful in deriving limits on decaying particles. Nevertheless, in the spirit of a conservative study we 
will also adopt 



^Li/^Li < 0.66 conservative. 



(6) 



as a second more conservative limit. Here this value derives by taking the average observed ^Li/^Li ratio of 0.042 and 
multiplying it by the ^Li/^Li destruction factor of 15.8. This is done with the understanding that any models which 
violate Eq. ([5]) but not Eq. ([6]) should be flagged as potentially being ruled out. On the other hand, such models, if 
not violating other observational constraints, could also explain the origin of the observed ^Li (cf. Ref. [SJ, |9|]). 

IV. CONSTRAINTS ON RELIC DECAYING PARTICLES 



In this section conservative constraints on the abundance of putative relic decaying neutral particles in the early 
Universe are presented. The number density nx of an early produced semi-stable particle species X with decay width 
Tx — 1/tx follows the equation 



dnx ,, „ nx 

at TX 



(7) 



10 




FIG. 9: Constraints on the abundance of decaying neutral particles as a function of life time for an Alx ~ 1 TeV particle with 
varying hadronic branching ratios. Constraints are shown for hadronic branching ratios (from bottom to top) log^^Bh ~ 0, 
—0.5, —1, -1.5, -2, -2.5, -3, -3.5, -4, -4.5, and -5, respectively, as labeled. The case Bh = is also shown. The labels always 
correspond to the solid (red) line above it. For each Bh the region above this solid (red) line is ruled out when conservative 
constraints are applied (i.e. Eas lH4l and[6|l. When the Li/ Li constraint Eq. [6]is replaced by the less conservative Eq. |5]the 
dotted (blue) constraint lines result. These dotted lines coincide for small and large tx with the solid (red) lines. 



^xh^ 




FIG. 10: As Fig. m but for Mx = 100 GeV. 
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where t is cosmic time and H is the Hubble constant. Here the first term on the right-hand-side describes dilution 
due to cosmic expansion, whereas the second term represents the particle decay. Results are presented for a large 
number of different hadronic branching ratios B^ of the particle X. In particular, it is assumed that a fraction Bfi of 
all decays are hadronic and a fraction 1 — B^ are electromagnetic. During hadronic decays the primaries injected into 
the plasma are assumed to be a quark-antiquark pair of total energy Mx and total momentum zero in the cosmic 
rest frame. The hadronization of the q-q fluxtube is followed with aid of the code PYTHIA, resulting in numerous 
nuclcons, mesons, e and 7-rays injected into the plasma. Note that for large Mx, as considered here, the spectrum 
of the resulting post-hadronization particles is virtually independent of the type of initial quark injected. As already 
stated in Sec. II, for sufficiently massive particles X the effects of electromagnetic decays on BBN yields depends 
only on the total energy Ee^j in e^ and 7-rays injected into the plasma. For electromagnetic decays Ef^ = Mx 
is thus assumed with no injection of nucleons or mesons. Hadronic decays are also associated with the release of 
energetic electromagnetically interacting particles. Following results from PYTHIA E^ — 0A5Mx is taken. Finally, 
it is noted that only baryon number conserving decays are considered. 

In Figs OS] constraints on the abundance of Mx = 1 TeV relic decaying particles for hadronic branching ratios 
Bh = 1, Bh = 3.333 X 10^^, and B^ = 10^'^, respectively, are shown. The colored (shaded) regions are excluded, 
with different colors (shades) corresponding to constraints derived from different light elements, as indicated in the 
figure captions. It is evident that the most stringent constraints derive from an overproduction of ^He at early times 
Tx ^ lO^sec, an overproduction of ^H in the decay time interval lO^sec ^ tx ^ lO'^sec, an overproduction of ^Li 
in the decay time interval lO^sec ^ tx ^ lO^sec, and an overproduction of the ■^He/^H-ratio for large decay times 
Tx ^ lO'^sec. The yellow area in the figures corresponds to elevated ^Li production in the range 0.1 < ^Li/'^Li < 0.66 
and should currently not be considered ruled out. Though in excess of the observationally inferred ^Li/^Li abundance 
ratio in low-metallicity stars it is conceivable that some of the produced ^Li has been destroyed in low-metallicity 
stars. The reader is referred to Sec. III.E for details. In contrast, as the source of the unexpected large abundance of 
^Li in low-metallicity stars is unknown, the yellow region may be considered as cosmologically interesting [8|, l9|. 

In Figs [Ql fTOl constraints on the abundance of relic decaying particles for a large number of hadronic branching ratios 
and two different particle masses, Mx = 100 GeV and AIx = 1 TeV are shown. The region above the solid (red) 
lines for each particular Bh are excluded. The figure also shows by the dotted (blue) lines how the constraints would 
change if one were to replace the observationally inferred limit on ^Li/^Li in Eq. [6]by the less conservative constraint 
Eq. [5] However, as argued above, the latter constraint should currently not be applied. 

V. CONCLUSIONS 

BBN is a powerful probe of the conditions in the early Universe and may severely limit the putative existence of 
relic decaying particles. In this paper, the effects of electromagnetically and hadronically decaying particles on the 
BBN light-element synthesis of ^H, ■^He, *He, ^Li, and ^Li have been considered in detail. All processes which the 
author believes to play a role at the 10% level have been included. These include nucleon-nucleon elastic and inelastic 
scatterings, nucleon-antinucleon annihilation, spallation of ^He by energetic nucleons, fusion reactions of ■^H and ■^He 
on He nuclei, mesonic charge exchange. Coulomb- and magnetic- moment scattering, Thomson scattering, Bethe- 
Heitler pair production, photon-photon scattering, (inverse) Compton scattering, and photodisintegration of light 
nuclei, among others. The analysis attempts to stay as close as possible to existing experimental data. A number of 
expressions, as for example the Coulomb stopping power of comparatively slow nuclei have been newly evaluated. All 
these non-thermal reactions are coupled to the thermonuclear fusion reactions known to be of importance during BBN. 
Results of such complete calculations have been first presented in Ref. [8,^] and particularly, Ref. ^], nevertheless, 
only for a limited number of decaying particle properties. 

Predicted light-element yields are compared to observationally inferred limits on primordial light-element abun- 
dances to infer accurate and conservative limits on electromagnetically and hadronically decaying particles in the 
early Universe. Constraints on the abundance of such decaying particles are shown for a large number of hadronic 
branching ratios and particles masses. This was necessary to obtain comparatively accurate limits as a simple extrap- 
olation between the results of only two hadronic branching ratios may yield to erroneous or inaccurate results. The 
results in this paper may be used to further constrain the evolution of the very early Universe. 



TABLE I: 
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No. 


Reaction 


Channel 


Ref. 


No. 


Reaction 


Channel 


Ref. 


1 


p + e* ^ elastic 


Coulomb 


App. D.l.a 


48 


p + p^ 


clastic 


[45, 46], App.E.3 


2 


p (N) + 7CMBR ^ 


Thomson 


App. D.l.b 


49 


p + "He ^ 


elastic 


[45], App.E.3 


3 


n + e='= -^ 


magnetic moment 


App. D.l.c 


50 


n + p ^ 


elastic 


[45,46], App.E.3 


4 


7 + 7CMBH, -^ 


e" +e+ 


App. D.2 


51 


n + "He-. 


elastic 


[45], App.E.3 


5 


e + 7CMBR -^ 


inverse Compton 


App. D.2, D.2.b 


52 


p + p^ 


p + p + 7r° 


[47], App. E.4 


6 


7 + p(^He)^ 


p(''He) + e" + e+ 


App. D.2. a 


53 




n + p + TT 


[47, 48], App. E.4 


7 


7 + e" ^ 


Compton 


App. D.2.C 


54 




2H + 7r+ 


[47, 48], App. E.4 


8 


7 + 7CMBH, — > 


7 + 7 


App. D.2.d 


55 
56 




P + P + tt+tt" 
^H + ^+^o 


[47], App. E.4 
[47], App. E.4 


9 


7 + 'H^ 


p + n 


[33 


57 




2H + 27r+7r- 


[47], App. E.4 


10 


7 + 'H^ 


m + n 


[37] 


58 




p + p + n'^TV~n 


[47], App. E.4 


11 




2n+p 


[37] 


59 




n + p + 27r+7r" 


[47], App. E.4 


12 


7 + ''He ^ 


^B+p 


[37] 


60 




p + p + 27r" + (7r"s) 


[47, 48], App. E.4 


13 




2p + n 


[37] 


61 




n + n + 2-K+ + (7r"s) 


[47], App. E.4 


14 


7 + ''He ^ 


^H+p 


[37] 


62 




n + p + 7r+ + (7r°s) 


[47, 48], App. E.4 


15 




^He + n 


[37,38] 


63 


n + p ^ 


p + n + 7r" + (7r°s) 


[47], App. E.4 


16 




^H + ^H 


[37] 


64 




2H + ^" + (^'M 


[47], App. E.4 


17 




^H + Tl+P 


[37] 


65 




P + P + tt" 


[47], App. E.4 


18 


7 + «Li ^ 


*He + n + p 


[37,39] 


66 




n + n + TX^ 


[47], App. E.4 


19 




^A + X 


[37] 


67 




n + p + n~n'^ 


[47, 48], App. E.4 


20 


7 + ^Li ^ 


*He + 3h 


[37,40] 


68 




p + p + n-n° + {n°s) 


[47], App. E.4 


21 




«Li + n 


[37] 


69 




^}i + 7V-n+ 


[47], App. E.4 


22 




■*He + 2n + p 


[37] 


70 




n + n + n+ + {n°s) 


[47, 48], App. E.4 


23 




*He + ^H + n 


App. D.3 


71 




n + p + 2n'2n+ 


[47|, App. E.4 


24 




"He + p^^Li+p 


App. D.3 


72 




n + p + n'n'^Tv'^ 


[47], App. E.4 


25 




23H+P 


App. D.3 


73 




p + p + 2n-n+ 


[47], App. E.4 


26 




^H + ^He + n 


App. D.3 










27 


7 + ^Be ^ 


"He + ^He 


[37,44] 


74 


p + "He ^ 


■■^H + 2p + {ns) 


App. E.5,E.6 


28 




^Li + p 


[33 


75 




^He + n+p+ (tts) 


App. E.5,E.6 


29 




*He + 2p + n 


[33 


76 




^H + 2p + n + (tts) 


App. E.5 


30 




"He + ^H+p 


App. D.3 


77 




^He + ^H 


App. E.5,E.6 


31 




^Be + n ^ "He + 2p + n 


App. D.3 


78 




""H + ^H+p+ins) 


App. E.5 


32 




23He + n 


App. D.3 


79 


n + "He -. 


^He + 2n + (tts) 


App. E.5,E.6 


33 




^H + ^He + p 


App. D.3 


80 
81 




^H + n+p+ (tts) 
^H + 271 + p + (tts) 


App. E.5,E.6 
App. E.5 


34 


7r+ +n ^ 


7r0(7)+p 


[15], App. E.l 


82 




^H + ^H 


App. E.5,E.6 


35 


tt" +P -^ 


7r°(7) + n 


[15], App 


E.l 


83 




^H + ^H + n + (tts) 


App E.5 


36 


R- +n-* 


p + X 


[15], App 


E.l 


84 


p + ^Li^ 


^He + "He 


App. E.5 


37 




n + X 


[15], App 


E.l 










38 


K- +p^ 


n + X 


[15], App 


E.l 


85 


^H + "He ^ 


®Li + n 


App. E.6 


39 




p + X 


[15], App 


E.l 


86 


^He + "He ^ 


'^Li + p 


App. E.6 


40 


Kl + u^ 


P + X 


[15], App 


E.l 


87 


"He + "He -^ 


•^Li + 2H ("Li +p + n) 


App. E.6 


41 




n + X 


[15], App 


E.l 


88 




^Li + p (■'Be + n) 


App. E.6 


42 


Kl+P—^ 


n + X 


[15], App 


E.l 










43 




P + X 


[15], App 


E.l 


89 


^He + p^ 


elastic 


App. E.6 










90 


'H+p-> 


^He + n 


App.E.6 


44 


p + p^ 


tt's 


[15], App. E.2 


91 


^H+p^ 


elastic 


App. E.6 


45 


n + p -^ 


tt's 


[15], App. E.2 


92 


^H^ 


^He + e- +17^ 


\49}, App.E.6 


46 


p + n ^ 


tt's 


[15], App. E.2 


93 


^Be + e" -^ 


^Li + Ue 


[49] 


47 


n + n ~* 


tt's 


fl5], App. E.2 
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APPENDIX A: GENERAL CONSIDERATIONS 

For the computation of light-element abundance yields in the presence of decaying or annihilating relic particles 
during or after BBN one has to essentially treat the thermalization process of energetic nucleons, antinucleons, 
and 7-rays. These latter particles result as primaries from the decay/annihilation process. Their thermalization 
occurs due to a host of processes discussed in detail in App. D and E. At the same time, thermal nuclear reactions 
of the thermal baryonic component may play a role at high temperatures T ^ 10 keV. In general, the problem 
simplifies considerably due to two facts: (a) thermalization may be treated as an instantaneous process on the Hubble 
time scale and (b) interactions between two non-thermal nucleons/nuclei are improbable. Given a typical proton 
density Up « 1.6 x 10^^/cm^(r/10keV)'^ and a typical nucleon ( nuclei )-nucleon interaction cross section of 30 mb 
and 300 mb at nucleon (nuclei) kinetic energy _Bfc « 1 GeV and 10 MeV, respectively, one finds a typical scattering 
time tn ~ 0.07sec(r/10keV)~^ which is approximately constant with kinetic energy. This may be compared to 
the Hubble time th « 1.4 x 10'*(r/10keV)-2sec to yield a ratio 10tn/th w 5 x 10-^{T/10keV)-^ where it is 
approximated that it takes around ten scatterings for the complete thermalization of a nucleon. It is seen that only at 
very low temperatures T < 0.1 eV thermalization may not be regarded instantaneous on the Hubble scale. Due to the 
short survival times of energetic nucleons (nuclei) only for excessively large ratios of non-thermalized-to-thermalised 
nucleons rinon-thcrmai/^6 ^ 10^(T/10keV) interactions between two non-thermalised nucleons become important. 
For such large f^non-thcrmai, however, BBN (almost) always fails badly in predicting the observationally inferred 
primordial abundances, such that such interactions may be neglected. 

Abundance yields due to decaying particles may be often understood as the product of individual factors. For 
example, for the production of ^Li nuclei by the fusion reactions 85 and 86, induced by non-thermal mass-3 nuclei, 
one may schematically write the final number of produced ^Li nuclei ANei^i, produced during a certain cosmic 
expansion interval d In a 



' d Decays \ / np,^ \ / ^3 \ fne 
d InR J \ Decay / V "•p,« / V "3 



ANsu = dNeu = dlna :i^f^ -^ -I:^ - 1 (Al) 



where a denotes the cosmic scale factor. Here the integrand is composed of the number of relic particle decays per 
logarithmic scale factor interval, times the number of energetic nucleons produced per decay, times the number of 
energetic '^H and '^He produced during the thermalization of a nucleon, times the fraction of '^H and '^He nuclei which 
fuse to form ^Li. Similarly there are contributions from photodisintegration not shown in Eg. lAll The present analysis 
treats in detail the third and fourth factors, pertaining to thermalization. For example, the last factor in the integrand 
(ne/ns) may be seen in Fig. [2] for various cosmic temperatures. The second factor is given by an assumed initial post 
decay state (e.g. qq an energetic quark- antiquark fiuxtube at an assumed energy) in combination with the results of 
an hadronization code, such as PYTHIA, which computes the number and energy of produced (anti)nucleons, mesons, 
and 7-rays after hadronization of the initial qq-fiuxtube. This factor (np^„/Decay) convolved with (n^/np^n) may be 
seen in Fig.[3](cf. also to Fig. |4] for a similar convolution due to photodisintegration). Finally, the first factor is given 
by the assumed abundance and life time of the putative relic. 

All figures in the main text limit Q-xh"^, the contribution of the decaying particle density to the present critical 
density fix if the particle were not decaying. Here h is the present Hubble constant in units of 100 km s~^Mpc~^. 
Other studies often present limits on the number of X-particles per photon times Mx , the mass of the X-particle, or 
the number of X-particles per radiation entropy times Mx ■ The conversion factors between these quantities are given 
by ^1 = nxMx/n^ = 2.5784 x 10"* GeVVixh'^ and ^2 = nxMx/s = 3.6639 x 10"^ GeYnxh"^- 

APPENDIX B: NUMERICAL MONTE-CARLO PROCEDURE 

The present analysis utilizes a Monte-Carlo approach to the problem. At each time step of the order of 10'* non- 
thermal nucleons and 6 x 10^ 7-rays are injected into the thermal plasma. Here there initial energies are randomly 
drawn, in the case of nucleons, from a kinetic energy distribution function generated by PYTHIA, and in the case of 
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photons, from the distribution function given in Eq. D.8. Mesons, which to a large part thermahze before interacting 
are injected in their number ratios predicted by PYTHIA. The subsequent evolution of each nucleon is followed by a 
Monte-Carlo (random probabilistic) sequence of events. Energy losses of nucleons are determined by a competition 
between nuclear scattering (spallation) processes (App. E.3 - E.5), in which a nucleon may loose a good part of 
its energy, and continuous energy losses, such as Coulomb stopping, Thomson scattering, and magnetic moment 
scattering (App. D.l). The probabihty of survival of a nucleon (nuclei) k against nuclear scattering while it looses 
continuously energy between energy Ei and Ef is given by 

where |d_E/d2;|c is the energy-dependent continuous energy loss per unit path length and 1/1% — 'Y^^cFki^...''n'i is the 
inverse of the nucleon (nuclei) mean free path against nuclear scattering. This probability distribution is mapped in 
the analysis by utilizing random numbers. Further, the evolution of any secondary energetic nucleons due to nuclear 
scattering is followed as well. The thermalization of nucleons is followed down to energies below which inelastic 
processes are not anymore possible. The remainder of the thermalization occurs quickly through elastic processes and 
does not further influence the BBN yields. In Fig. [T] an approximation to the integrand in Eq. IBll may be seen for 
protons and neutrons at various cosmic temperatures. 

Photons are treated in a similar way. Each injected photon and the by it produced secondaries are followed in 
a Monte-Carlo way until all photons have dropped below the lowest photodisintegration threshold. The number of 
nuclei (nucleons) A per decaying particle X produced due to all 7-rays created as a consequence of the decay of X 
are computed via 

D^ = (E ^^^0 (^ "^ ^ -B,->A...crAE^) - n, Y. -A,^c...crAE^)) ■ (B2) 

Here the sum over index i runs over the energy spectrum of the initial (after CMBR cascade) 7-ray spectrum, a are 
photodisintegration cross sections, c is the speed of light, and Tj{E-') are 7-ray survival times against the processes 
discussed in App. D.2.a, D.2.C, D.2.d, and D.3. The sums over i and j are different due to the inclusion of secondary 
photons. In particular, a photon with initial energy El^ may produce secondary photons j due to the processes 
described in App. D.2. The sum over j thus includes all 7's, primary and secondary, tertiary, ... as long as they 
are above the photodisintegration threshold. Individual survival times are generated by a Monte Carlo and the 7-ray 
survival probability P^{E~^) = exp{Tc/l^{E^)), with 1^{E^) the 7 mean free path. Finally, it is noted that for accuracy 
reasons the initial first-generation photon spectrum Eq. D.8 is sampled by an equal amount of photons above and 
below the photodisintegration threshold for ''He. This was necessary as otherwise much larger (^ 6 x 10'^) numbers 
of photons would have to be followed in order to obtain accurate result. 

APPENDIX C: KINEMATIC RELATIONS 

For the convenience of the author this appendix presents some kinematic relations required to convert cross section 
data presented in terms of center-of-mass (CM) quantities to cross section data in the laboratory (L) frame of reference. 
For the energy gain of particle 2 (initially at rest) in elastic 2^2 scattering, and given the scattering angle in the 
CM frame, one finds 

^El,2 = 7'cm«cmM2(1 - cosGcm) (CI) 

where vqm = p\/{M2 + Mi + E\^^^^, Ecm,2 = 7i,cm^2, and |pcm| = i^CMlfcM^s. (Note that |pcm| is invariant 
during elastic scattering.) Concerning inelastic endothermic 2^-2 scattering reactions 1-1-2^3-1-4 with 2 initially 
at rest (e.g. reactions 77, 82, 85-88) the threshold in the laboratory system is given by Ej^ = Am(l + 7711/1712) with 
Am = (777,3 + "^4 ~ "^1 ~ 1^2), the mass differences. During inelastic 2 — > 2 scattering with particle 2 initially at rest, 
\PCM I is not anymore invariant due to conversion of kinetic energy to rest energy. Utilizing the above given (before 
scattering) ucm , bcM | relations one may solve the implicit equation i^cM + -^cm — -^cm + ^cm for p^M 1 '^i*'^ -^'cm '^'^d 
^CM containing the after scattering modified p'^m- Provided one has 0cm the CM scattering angle, the laboratory 
energy of particles 3 and 4 are determined by i?£ — ^vcmi-Ecu + ''^CmIpcmI'^osOcm) and similarly for particle 4 with 

WCM -^ — WCM- 



15 



APPENDIX D: NON-THERMAL ELECTROMAGNETIC INTERACTIONS 

1. Electromagnetic Stopping of Non-thermal Nucleons 

a. Coulomb interactions 



Stopping of charged protons and nuclei by Coulomb scatterings off ambient electrons (and positrons) as well as 
excitation of collective plasma modes (involving a larger number of electrons) is particularly important. In the early 
Universe single particle scatterings off electrons and positrons dominate. Since many of these occur with only small 
momentum transfer (i.e. small-angle scatterings) the process may be regarded as a continuous energy loss for protons 
and nuclei. Following Jackson [50], Gould [51|, and own calculations the energy loss per unit path length traveled 
may be found 



dx 



-UJ„ 



where 



0.76v 



In 



47rn„± a 



1 



(Dl) 



(D2) 



is the square of the plasma frequency, and 



Zc 



1 



"fmeV^ jnieV 



(D3) 



In these expressions Z, v, and 7 refer to electric charge number, velocity, and Gamma-factor of the nuclei, respectively, 
?7ie is the electron mass, a the fine structure constant, and ne± the total number density of electrons and positrons. 
This latter is given by 



>l/2 



'e.nct 



e,pair> 



with 



1 / 2 



(!)"("-)"v--(-£E) 



(D4) 



where T is cosmic temperature. Here rig. not refers to the electrons needed for charge neutrality of the Universe. 
Eq. (jD4p is only accurate in the T/nie <C 1 limit. Note that the second term in Eq. (jDip derives from the magnetic 
moment scattering contribution of the electron absent in the Rutherford scattering cross section but present in the 
correct Mott scattering cross section. Note also that, in contrast to the findings of Ref. [15[ it is found that Eq. (jDip 
is also valid in the extremely relativistic limit up to 7 < Mp/rrie, where Mp is the proton mass. Nevertheless, 
an important modification to Eq. (|D1[) occurs for small nuclear velocities and large temperatures when the nuclear 
velocity v falls below the typical thermal electron velocity Ve- In this case the Coulomb stopping power decreases by a 
factor Ve/v^ 15|, a very larger factor 500 for, example, for tritium nuclei of 10 MeV kinetic energy at a temperature of 



T = 20 keV. In order to obtain an accurate result for the stopping length in this regime an adequate thermal average 
over the electron velocity distribution has to be performed, as Coulomb stopping is dominated by Coulomb collisions 
off slow electrons. After a lengthy calculation the author finds 



dEc _ ATT{Zaf 
dx 



rrieV 



n^± erf(y) 



V \ zm,. I 2 \ V mp 



4=e-^^fAfl--^ 



T 8 



T 



2m„ 



(P5) 



for the Coulomb stopping which includes a thermal average over a Maxwell-Boltzmann distribution of electron 
(positron) velocities. Here y = yJrripXp- I2T such that for y ^ 1 one recovers the result Eq. (|D1[) where the ve- 
locities of all electrons are below that of the nucleus. The term A in Eq. (|D5[) denoted the logarithm in Eq. (jDip . 
Note that in Eq. (jD5p an expansion has been performed to lowest non-trivial order in the small parameters v and 



y/r/rrif., where both have been regarded to be of the same order. 



Thomson scattering 



Energetic protons may loose energy to the CMBR by Thomson scattering, a process which is particularly efficient at 
high CMBR temperature. Due to the large number of CMBR photons many scatterings are involved and the energy 
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loss of protons due to Thomson drag may also be treated as a continuous energy loss. Following Reno & Seckel 15 1 
and own calculations one finds 



dE- 



Th 



da; 



~9 15~iW| 



2 r74 

IpVpZ 



(D6) 



for the energy loss per unit length where a and T are fine structure constant and CMBR temperature, Mp the proton 
mass, 7p and Vp its gamma factor and velocity, and Z its electric charge. Note that Eq. (|D6p is applicable for protons 
of arbitrary relativity. 



c. Magnetic moment scattering 

Neutrons may loose energy due to scattering off electrons and positrons by virtue of their anomalous magnetic 
moment. This process is only important at higher temperature T ~ 50keV due to the larger numbers of e . By aid 
of a calculation employing the relevant interaction amplitude one may find 



dg„ 
dx 






7n^lri^± : 



(D7) 



where n = —1.91 is the anomalous magnetic dipole moment and all other notation is as before. Eq. (jD7P applies for 
non-relativistic as well as relativistic neutrons. 



2. Electromagnetic Cascades Induced by 7-rays and Charged Particles 

The injection of energetic electromagnetically interacting particles in the early Universe leads to a cascade on the 
CMBR with 7-rays pair-producing (7 -I- 7cmbr -^ e+ -|- e~) and the produced electrons (positrons) inverse Compton 
scattering (e^-t-7cMBR — * e^-l-7) off the CMBR photons 52, l53|, |5J, |55|, |50| . This cascade is very rapid due to the large 
number of CMBR photons. Only when 7-ray energies have fallen below that for pair production E^ < ml/ Ecmbr 
interactions on protons and nuclei become important. The spectrum of these "breakout" photons below the pair 
production threshold depends essentially only on the total electromagnetically interacting energy Eq injected above 
Ec- It has been found to be well approximated by ISSllS 



diV^ 
dK 



' Ex 
' Ex 



-1.5 



for E^ < Ex 



for Ec>E^> Ex 



(D8) 



where Kq = Eo/{E'j^[2 + \n{Ec / Ex)]) is a normalization constant such that the total energy in 7-rays below Ec 
equals the total energy Eq injected. Following the analysis of Ref. [56| a value of Ec ~ rrv^/22T is employed in the 
present analysis with Ex ~ O.OSi^c [53]. These values are very close to those (Ec ~ ml/23.6T and Ex ~ 0.0264£'c) 
advocated by Ref. 55]. Subsequent interactions of these "break-out" photons include photon-photon scattering 
(7 + 7CMBR, ^7 + 7) mainly redistributing the energy of energetic 7-rays right below energy Ec, Bethe-Heitler 
pair production (7 -1- p(^He) — > p('*He) + e~ + e"*"), Compton scattering (7-1-6" — > 7 4- e~) off thermal electrons, 
with the produced energetic e~ inverse Compton scattering to generate further low-energy 7-rays, as well as nuclear 
photodisintegration. All these processes are included in the analysis and their detailed treatment is described below. 



a. Bethe-Heitler pair production 



The cross section for e pair production by 7-rays of energy E^ scattering off protons and helium nuclei at rest is 
given by 



3a /28, 

8 TT 



2E^ 



218X , 
27 / 



(D9) 
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where (iTh, a, 'me, and Z are Thomson cross section, fine structure constant, electron mass, and nuclear electric 
charge, respectively. The expression is valid in the regime 1 <C E^/nie <C a~^Z~^'^. For energies in the range 
2me < i?-y < 4 MeV the cross section is roughly approximated as constant at the value of Eq. (|D9[) a,t E^ — A MeV. 
The energy of the produced e~ (e+) is approximated by a flat probability distribution within the range [TTig, E^ — 277ie] 
with the antiparticle e"*" (e~) carrying the remainder of E-y. 

b. Inverse Compton scattering 

Electrons and protons produced during the dominant Bethe-Heitler process may upscatter CMBR photons to 7- 
ray energies by inverse Compton scattering. As these 7-rays may later photodistintegrate nuclei, a fairly accurate 
treatment of inverse Compton scattering is required. Here the exact inverse Compton scattering rate is unimportant 
since it is always large compared to the Hubble expansion. For the cumulative probability distribution that an e~ 
(e"*") of energy E^ = jerrie upscatters a CMBR photon to energy E^ = Aj'^Tx, where T is CMBR temperature, one 
may find 

P^c « e-^ (1 + J^ - \^') + r(0, x) Q.t3 - ^x) , (DIO) 



r(0,x)= / dye-y/y (Dll) 



with 



is an incomplete Gamma function. Note that for the evaluation of Eq. (jD10[) the approximation that the e* are 
ultrarelativistic, i.e. 7e ^ 1, and that the scattering process occurs in the Thomson regime, i.e. li^^T /me ^ 1, 
have been made. These approximations are appropriate since only ultrarelativistic e may produce 7-rays by inverse 
Compton scattering on the CMBR which are energetic enough to photodisintegrate nuclei, and since e* produced 
during Bethe-Heitler pair production are not energetic enough for the inverse Compton scattering to proceed in the 
Klein-Nishina limit. Furthermore for the derivation of Eq. (|D10p the Einstcin-Bosc occupation number {e~^'' — 1)~^ 
has been approximated by ^(3)e^^^ with i^(3) « 1.2021 such that the total number of CMBR photons remains 
the same. The approximation was required to find a closed form for Pic- By comparison to an accurate numerical 
evaluation one finds that the approximation is typically good to within 3% - 10%. Lastly one finds for the average 
scattered CMBR photon energy a value of (Ej) = Jei^CMBR) with (i^cMBR) ~ 2.701T the average CMBR photon 
energy. 

c. Compton scattering 
The cross section for 7 rays of energy E^ scattering off electrons at rest in the Klein-Nishina limit E^ ^ nie is 



given by [57 1 



3CTTh 



8 oj 



lii2uj + - + 0{ 

2 \ OJ / 



a;>l, (D12) 



where lo = E^Jm^ and (iTh is the Thomson cross section. The cumulative probability distribution that the 7 scatters 
to energy E' < E!r during the Compton scattering process may be derived by utilizing the results in Ref. (^57]) and 
is found as 

hx'^ - Ini + ln2w , , , , , 

PeiE^,E'^ < E^^) « ^ ,^'^^^^ + 0(l/c.) , (D13) 

where x = E'ijE^ and the kinematic limits of x fall in the range 1 and 1/(1 + 2lS) w 1/2w. The scattered electrons 
may produce further 7-rays by inverse Compton scattering on the CMBR, an effect which is taken into account in 
the calculations. 
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d. 7-7 scattering 



Following the analysis of Svensson & Zdziarski 53|, |54| the total rate for scattering of a 7-ray of energy E-^ on a 
blackbody of temperature T is given by 



R 



24l397r-^ 
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Employing Eq.(3.9) of Ref. [5j| one may compute the cumulative probability distribution for the energy of the scattered 
photon E' to be below E^ as 



10 



Pc{E^,E'<E')^—x 



1- x + x'^ x^ + -x'^ 

2 5 



(D15) 



where x = E'i^jE-^ < 1. The energy of the scattered blackbody photon is given approximately by E' bb ~ ^1 ^ ^'j 
in the limit T -^ E^ and its effects are treated in the calculations as well. 

3. Photodisintegration 

Photodisintegration of the light-elements by 7-rays generated as a consequence of relic particle decays plays a par- 
ticularly important role in the determination of BBN yields, as noted early on 1^. It is only effective at temperatures 
below T < 6keV (cf. Fig. [5]), since at higher temperatures 7-rays of sufficient energy Eg « 2MeV to photodisin- 
tegrate ^H, ^Li, and ^Be predominantly pair produce on the abundant CMBR photons (cf. App. D.2). Below this 
approximate cosmic temperature pair-production is kinematically forbidden such that photodisintgeration becomes 
probable. The competing Bethe-Heitler pair production, Compton scattering and 77 scattering were described in the 
previous subsections. 

The (lowest) photodisintegration thresholds for the various light nuclei synthesized during BBN 
2H,3H,3He,4He,6Li,^Li, and ^Be are given by 2.225 MeV, 6.257 MeV, 5.493 MeV, 19.814 MeV, 3.699 MeV, 2.467 MeV, 
and 1.587 MeV, respectively. The present analysis utilizes photodisintegration data as parametrized by Ref. [19| 
with a few of these parametrizations changed to produce improved fits to the available reaction data. For original 
references the reader is referred to the references in Ref. \1^. 

In addition, the reactions 23-26, and 30-33, not considered in Ref. [l9|, have been included in the present analysis. 
Here the following parametrizations of reaction rate data 42, ISa, ISSj have been adopted: 



Reaction 23 k2 



a{E^) = 3.8mbQ2-3(^^ _ Qa)/E^-'^ 



2.1mhQl-^{E^ - Qi)/E^-^)e{E ~ Qi) 



(D16) 



with Qo = 8.725 MeV and Qi = 23 MeV and Q{x) the step function. 
Reaction 24 ^\E^: 



(7{E^) = 10.8mbQ2(^^ _ Qy-^/E^ 



(D17) 



with Q = 9.98 MeV, and 
Reaction 25 [H: 



a{E-y) = 1.44 X l0^mhQ^"{E-y - Qf-^/E^ 



(D18) 



with Q — 22.28 MeV. Reaction 26 has been parametrized as reaction 25 but with Q ~ 23.05 MeV since no experimental 
data exists. Similarly, since no reaction rate data exists for the cross sections of reactions 30-33 they have been 
approximated by the data of their respective mirror reactions 23-26 but with binding energies appropriately replaced 
by, Q = 7.08 MeV, Q = 10.68 MeV, Q = 22.17 MeV, and Q = 21.40 MeV for reactions 30,31,32, and 33, respectively. 
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APPENDIX E: NON-THERMAL HADRONIC INTERACTIONS 
1. Mesonic Induced Charge-exchange reactions 

Injection of mesons in the form of metastable pions and kaons may induce the charge- exchange reactions 34-43 
to convert protons to neutrons and vice versa [l^. As essentially each neutron is incorporated into a helium nuclei 
at r « 80 keV, mesons may thus effect the ^He abundance. Since more protons are converted to neutrons than 
vice versa, injection of mesons leads to an increase in the ^He abundance. Mesonic charge-exchange reactions are 
particularly important at fairly high temperatures. At T ~ 1 MeV, for example, a fraction ~ 0.01 of all mesons induce 
charge exchange with the remainder decaying. In contrast, this fraction drops to ~ 10^^ at T « 100 kcV. The present 
analysis follows the treatment of Ref. 'l5-| and the reader is referred to this study for details. 

The mesons may be due to a variety of processes. It is well known that the hadronic decay (or annihilation) results 
in a multitude of mesons. However, under certain circumstances even electromagnetic- and weak- decays may yield 
mesons. This may be due to decays involving r^ in the final state, which have a branching ratio to pions. Similarly 
injection of high-energy E^, ^ 100 GeV neutrinos at T '--^ IMeV may pair-produce tt on the neutrino background. 



Pair-production of pions by 7-scattering on the CMBR [2l| is unlikely as the Compton cross section is typically a 
factor ^ 300 larger. 

Production of mesons by elcctromagnetically or weakly interacting particles has not been included in the present 
study. In the case of hadronic decays their effect is subdominant, and in the case of electromagnetic- and weak- decays 
the analysis is model-dependent on the decay products. 

2. Nucleon-Antinucleon Annihilations 

Hadronic decays lead to the production of nucleons and anti-nucleons. It is assumed here that the total baryon 
number is conserved during the decay. At higher temperatures (T > 90keV for n and n and T > 20keV for p and 
p) they thermalize by electromagnetic energy losses (cf. App. D.l) with antinucleons annihilating at rest thereafter. 
The combination of injection of nucleons (with a predefined n/p-raiio, typically close to one) and annihilation of 
antinucleons on the pre-existing nucleons may lead to an increase of the neutron-to-proton ratio and the '*He abundance 
resulting from BBN [l5|. Antinucleons mostly annihilate on protons due to the small thermal n/p-iatio and due to 
Coulomb enhancement of the pp annihilation cross section. The annihilated protons are replaced by the injected 
nucleons which typically come in a ratio n/p ~ 1 thus yielding an effective increase of the n/p-raiio. Given a typical 
injected meson-to-(anti-)baryon ratio of ~ 20 in hadronic decays the effects of nucleon anti-nucleon injection are far 
more important than those of mesonic charge exchange reactions. This holds particularly true at lower temperatures 
T < 300 keV. The present study uses the reaction rate data as given in Ref. [I3 . 

When electromagnetic stopping becomes less dominant at lower temperatures thermalization of energetic nucleons 
occurs partially due to nucleon-nucleon scattering and nuclear spallation (cf. App. E.3,E.4,E.5). Antinucleons, on 
the other hand, never thermalize but rather annihilate during one of their first interactions with thermal protons or 
helium- nuclei. As an energetic nucleon has of the order of ^ 20 N — N scattering events before dropping below the ^He 
spallation threshold it's probability to inflict a '''He spallation event is factor ^ 10 larger than that of an antinucleon. 
It has been therefore approximated that antinucleons do never spall ^He, leading to an about 10% underestimate in 
^H, '^H, and '^He production due to particle decay. 

3. Elastic A'^-A*' scattering 

Elastic nucleon-nucleon scattering and nucleon-^He scattering of energetic neutrons and protons off the thermal 
protons and ^He is important in the nucleon thermalization process in particular at lower temperature [T < 90keV 
for n's and T < 20 keV for p's) and when the nucleons are not too energetic, i.e. E ^ 1 GeV. This may be seen 
in Fig. [TJ In contrast, very energetic nucleons, E ^ 1 GeV, are stopped at lower temperatures mostly by inelastic 
scatterings of protons and ^He (cf. App. E.4, and E.5) and at higher temperatures either by Thomson scattering off 
CMBR photons in the case of protons (cf. App. D.l.b) or magnetic moment scatterings off thermal e^ in the case of 
neutrons (cf. App. D.l.c). Elastic scatterings off protons also produce secondary energetic protons which are taken 
into account in the present study. 

For the elastic p-p, n-p, p-^He, and n-*He scattering cross sections we take the data as compiled by Ref. |45|. This 
compilation compares well with that given in the more recent compilation given in Ref. 49]. It is noted here that 
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due to approximate isospin invariance the p-^He and n-'*He cross sections are almost equal. To calculate the energy 
loss in p-p and n-p scattering events, an against nuclear data tested algorithm presented in Ref. [46| is used. Each 
scattering event is treated by a Monte Carlo. This algorithm predicts an almost uniform probability distribution for 
fractional energy loss AEk/Ef- for < 250 MeV protons progressing smoothly towards a bimodal probability distribution 
(either forward scattering, i.e. AEk/Ek « or backward scattering, i.e. AEk/E^ ~ 1) when the p kinetic energy 
is increased towards 2 GeV. Neutron-proton scattering follows a similar though less well-defined trend. For kinetic 
energies E^ > 2 GeV the probability distributions of 2 GeV nucleons are utilized. 

Nucleon-^He scattering follows a very different pattern with '*He recoil energies to high probability very small. 
Given data by Ref. [60|] the typical ^He recoil energy is approximated at 5 MeV in each event. It has been verified 
that a change in this recoil energy has negligible effect on the BBN yields. For a more detailed discussion on the ^He 
recoil energies which play a role in ''He-^He fusion reactions 87 and 88 the reader is referred to App. E.6. 

4. Inelastic A'^-A'^ scattering 

For nucleon-nucleon scattering with nucleon kinetic energies in the GeV range scattering processes are mostly 
inelastic accompanied by the production of pions. Treating these processes properly is important not only for a 
determination of the typical nucleon energy loss, but also because scattering processes may easily convert protons to 
neutrons and vice versa. In fact, due to a large asymmetry in the cross sections the conversion of energetic protons to 
neutrons occurs more often than the inverse of this process. This has implication for spallation processes, increasing 
^H,'^H,^He, and ®Li yields, as neutrons are stopped by scattering off protons and ^He for temperatures T ^ 90 keV, 
whereas Ep ^ 1 GeV protons mostly loose their energy via Coulomb stopping. This study uses the inelastic N-N 
cross sections reactions 52 - 73 in Table 1 as compiled by Ref. |47| , where some of the proposed fitting functions had 
to be modified by the author the correctly account for the data. Concerning the total average energy loss /^r into pion 
mass and kinetic energy in inelastic scatterings experimental data by Ref. [6l| at Ep = 790 MeV and a simulation 
with Pythia at Ep — 53 GeV has been fitted logarithmically to obtain /^ « 0.2 -|- 0.081n(i?/0.8GeV) The remainder 
of the energy has been split in proportions (0.15,0.85) between the outgoing nucleons roughly equal to the average 
energy split for energetic elastic collisions. For the nucleon energy loss in the inelastic spallation processes discussed in 
App. E.5 it has been assumed that the outgoing nucleons are missing the binding energy necessary for the spallation 
process, and the remaining energy is shared between primary and secondary nucleon in proportion (0.75,0.25). This 
corresponds to the average values for elastic N-N scattering in the several hundred of MeV range. 

5. Nuclear spallation 

Aside from ''He photodisintegration which is operative only at lower cosmic temperatures T < 3 keV '''He spallation 
reactions at temperatures T < 90 keV induced by energetic nucleons are the main source of additional ^H, '^H, ■^He, as 
well as ^Li (through nonthermal fusion reactions induced by energetic 3-nuclei, cf. App. E.6) in BBN with decaying 
particles. It is therefore important to employ accurate reaction cross section data. The ''He-spallation reactions 
74 - 83 in Table 1 have been included in the analysis. Experimental data on these reactions has been obtained at 
nucleon energies E = 28 MeV [^ , 53 MeV , 90 MeV ^ , 220 MeV ^ , 300 MeV [^ , 620 MeV [67] , 1 .42 GeV \E^ , 
and 2.61 GeV |68| , respectively, and has been interpolated in between. This data and the interpolation is shown in 
Fig. [TlJ The data also includes the possible production of additional pions in the spallation process. Note that 
^He-spallation data induced by either neutrons or protons has been included on equal basis in the compilation. Here 
reactions induced by neutrons, such as *He(n,np)'^H, have been associated with their corresponding mirror reactions 
induced by protons, i.e. ^He(p, np)'^He. This is theoretically and experimentally justified due to approximate isospin 
invariance, as long as Coulomb effects are negligible, in particular, for nucleon energies larger than a few MeV. For 
the energies of nucleons and nuclei's after the spallation the reader is referred to App. E.4 and E.6. The spallation of 
energetic freshly synthesized ^Li reaction 84 has also been included [69| but plays only a minor role for the results. 

6. Non-Thermal Fusion Reactions 

Non-thermal fusion reactions are particularly important for the formation of ^Li which during thermal BBN is 
produced only in very small amounts. When relic particles decay during or after BBN ^Li is therefore the first element 
which becomes significantly perturbed compared to its SBBN abundance. This has been shown first in the case of 
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FIG. 11; Compilation of the available *He spallation data as a function of nucleon kinetic energy (with ''He at rest) for the 
spallation reactions 74,75,76, and 78. The lines joining independent data points (74 - solid; 75 - long-dashed, 76 - short-dashed, 
78 - dotted) represent the fit used in the present analysis. Refer to the text for references. 



hadronic decays in Ref. [lj| and later in the case of electromagnetic decays in Ref. 17]. The most efficient reaction 
sequence is initiated by either the spallation *He(A^, ... or photodisintegration ''He(7, ... of *He to produce energetic ^H 
and ^He nuclei. These later may then fuse on ambient ^He-nuclei to form ^Li via •^H(^He, n)^Li and •^He(''He,p)^Li. If 
only a small fraction ~ 10~'* (cf. Fig. ^ of these energetic 3-nuclei fuse an observationally important ^Li abundance 
may result. Energetic 3-nuclei in the tens of MeV range are essentially all stopped by Coulomb interactions. Since this 
process is less efficient at higher energy (i.e. (1 / E-i)dE-i / dx oc E^^] cf. App. D.l.a) and since the fusion cross section 
does not drastically fall with energy the formation of ®Li is very sensitive towards the high-energy tail of produced 
3-nuclei. It will be seen that the lack of knowledge of this tail actually introduces some uncertainty in the ^Li yield. 

In this paper the ^H('*He, n)^Li and '^IIe("'He,p)^Li cross sections as determined with the aid of reverse reaction 
rate data by Ref. [70| and Ref. [19|, respectively, have been adopted. Here the first reaction has a threshold in the ''He 
rest frame of Eth — 8.39 MeV, whereas the second reaction has Eth — 7.05 MeV. Cross sections for the spallation and 
photodisintegration reactions are taken as described in App. E.5 and App. D.3, respectively. Concerning the "recoil" 
energy distributions of the 3-nuclei in the case of photodisintegration it is given by E^ = {E^ — E J^'f^)'mff / {mp^ +7113) for 
•^H nuclei, whereas for spallation one has to resort to experimental data. Data has been published for the ^He(n, np)^H 
reaction 1641 at En = 90 MeV in the ''He rest frame and more recently for the ^He(p,pp)'^H and ''He(p, np)'^He 
reactions |7lj at Ep — 220 MeV. Data at higher energies exists as well 72| which lies above the threshold of pion 
production during the spallation process. In general, the energy distributions of 3-nuclei seem almost independent 
of incident nucleon energy though care has to be taken in the comparison since sometimes only partial (quasi-free- 
scattering) distributions are shown (as in the case of Ref. [72|). In this study the complete distributions as given 
in Ref. 71| are utilized (in binned form) for both the '*He(p, ... reactions and their '*He(n, ... corresponding mirror 
reactions. The utilization of reaction rate data of the corresponding mirror reactions seems not only theoretically 
justified 45|, |73| but also when experimental data of Ref. 64] and Ref. [7l| are directly compared. 

Some uncertainty in the ^Li yield results from an inaccurate determination of the few percent high-energy tail in the 
3-nuclei ener gy d istribution function. The highest energy events are particularly difficult to detect in bubble chamber 
experiments |73| , such that the distribution of Ref. [71| only extends towards E3 w 90 MeV. In Ref. ^M] a fraction 
of about 2% of all spallation reactions seem to result in '^H nuclei at the kincmatically highest allowed energy. This 
corresponds to backward scattering and such peaks in the distribution functions arc also observed in '*He-'*He elastic 
scatterings [74| . If the four-event backward scattering peak as observed in Ref. [64] is real, rather than a fluctuation 
or measurement error, resulting ^Li yields may increase by some tens of per cents due to these high-energy 3-nuclei 
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(cf. Fig. [2]). However, on conservative grounds this study neglects ^Li production due to the very highest energy 
3-nuclei. 

Energetic 3-nuclci may also be produced by *He(n, ^H)'^H and its mirror reaction. For the e nerg y distributions of 



this subdominant 3-nuclei channel reaction data of the compilation/determination in Refs. 45|, |75| at proton kinetic 
energies of 55 MeV, 85 MeV, and 156 MeV, respectively, have been interpolated. Fusion reactions such as ^H(a, 7)®Li, 
'^H(a,7)^Li, and ^He(a,7)'^Be have not been considered as cross sections involving the emission of a photon are 
typically only in the nb to ^b range. The N{N, 7r's)^H reaction, on the other hand is treated as discussed in App.E.4. 

Note that in a narrow temperature interval around T ~ 3 — 20 keV Coulomb stopping of not too energetic E3 ~ 
10 — 30 MeV 3-nuclei looses its efficiency (cf . App. D.l.a) due to the nuclei velocity falling below a typical electron 
thermal velocity. In this regime of efficient ^Li synthesis other processes reactions 89 - 91 which contribute to the 
thermalization of 3-nuclei may become of some limited importance. They have been included in the calculation by 
taking data from Ref. [76[ . At very low temperatures T < 1 eV energetic '^H stopping times may become longer than 
the "^H decay time. It is therefore possible that '^H decays to ^He before it falls below the ^Li fusion threshold and 
this effect has been taken into account. 

Fusion reaction between energetic ''He (produced by iV-^'He elastic scatterings) and thermal '^He-nuclei may, in 
principle, also lead to the synthesis of ^Li, as well as ^Li. In practice, however, this process is subdominant and has 
therefore not been included. The threshold energies for ^He(a,p)^Li and '^}ie{a,pn)^Li are 37.8 MeV and 49.2 MeV, 
respectively. These are much larger than those for the fusion of 3-nuclei on ^He. Given the data in Ref. [60| 
the distribution function for recoil energies of '*He-nuclei in A^-^He elastic scattering may be approximated by an 
exponential distribution dP/dEn « aexp(— 6£^4) with a « 225mb/GeV and b k, 65GeV~^, independent of nucleon 
kinetic energy. The total cross section for producing E/^ > 50 MeV recoil ''He in A-^He scatterings is thus in the 
^ 0.1 mb range. Note that E^ >50MeV has been chosen as only for such energies the ^He-^He fusion cross section 
becomes appreciable ~ 100 mb [77|. The cross section of ~ 0.1 mb for the production of ''He nuclei energetic enough 
to produce ^Li by fusion has to be compared to the typical total cross sections for ^He(A^, ...)'^H, etc. of ~ 30 mb times 
the fraction ^ 0.3 of 3-nuclei energetic enough to fuse to ^Li. This yields '--' 10 mb compared to 0.1 mb. Nevertheless, 
^Li fusion via ^He-nuclei is enhanced due to a less efficient Coulomb stopping at 50 MeV versus the canonical 10 MeV 
for 3-nuclei, resulting in an enhancement factor of (50/10)^/4 « 6. Here the 4 in the denominator comes from the 
factor two higher charge of ''He nuclei compared to ^H nuclei. Altogether it is found that energetic ^He contributes 
only of the order ~ 6% to the final ^Li abundance This is also not changed by the existence of a sharp backward 
scattering peak in A^-'^He scattering [7^ which, similarly to the possible existence of a backward peak in the 3-nuclei 
distribution (see above) may increase the ^Li synthesis due to energetic ^He by a factor ^ 1.2. 

APPENDIX F: THERMAL NUCLEAR REACTIONS 



Thermal nuclear reactions are treated by the Kawano code updated by the NACRE [4]J reaction compilation. For 
the neutron life time the value t„ — 885.7 sec is used as quoted as the 'world' average in the Particle Data Booklet. For 
the fractional contribution of baryons to the critical density the central value as determined by three years of WMAP 
observations [4] fibft,^ = 0.02233^Qgg^ is applied. Given the conversion Vlbh^ = 3.65 x 10^77 this corresponds to a 
baryon-to-photon ratio of 77 = rib/n^ — 6.12 x 10"""'. With these input values the abundances within a SBBN scenario 
are found at ^H/H = 2.68 x 10-^ ^He/^H = 0.38 (^He/H = 1.04 x 10~^) Yp = 0.248, and ^Li/H = 4.34 x lO^i^ 
respectively. The ^Li abundance in a SBBN scenario is negligible compared to the observationally inferred value of 
^Li in low-metallicity stars. 
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